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In previous work [Amendola et ai, Phys. Rev. D86 (2012) 063515], Principal Component Analysis 
based methods to constrain the dark energy equation of state using Type la supernovae and other 
low redshift probes were extended to spectroscopic tests of the stability fundamental couplings, 
which can probe higher redshifts. Here we use them to quantify the gains in sensitivity obtained 
by combining spectroscopic measurements expected from ESPRESSO at the VLT and the high- 
resolution ultra-stable spectrograph for the E-ELT (known as ELT-HIRES) with future supernova 
surveys. In addition to simulated low and intermediate redshift supernova surveys, we assess the 
dark energy impact of high-redshift supernovas detected by JWST and characterized by the E-ELT 
or TMT. Our results show that a detailed characterization of the dark energy properties beyond the 
acceleration phase (i.e., deep in the matter era) is viable, and may reach as deep as redshift 4. 


I. INTRODUCTION 

The observational evidence for the acceleration of the 
universe demonstrates that our canonical theories of cos¬ 
mology and fundamental physics are at least incomplete, 
and possibly incorrect. Substantial efforts are being put, 
therefore, into identifying the mechanisms responsible for 
it—see for a recent review. At the most basic level 

one would like to know if the acceleration is due to a cos¬ 
mological constant (the simplest viable possibility), to a 
modification of the behavior of gravity on large scales, or 
to a cosmological dynamical scalar field, but ultimately 
the goal must be to characterize its behavior as a func¬ 
tion of redshift. Mapping the redshift dependence of the 
dark energy equation of state is a simple, standard way 
to do this. 

Astrophysical tests of the stability of fundamental cou¬ 
plings, such as the fine-structure constant a, are among 
the most powerful probes of the equation of state of dy¬ 
namical dark energy Si- The idea was qualitatively 
discussed in [3-[l], and more recently forecasts of its im¬ 
pact have been obtained 1,0- This complements other 
methods due to its large redshift lever arm and the fact 
that these measurements can be done from ground-based 
facilities, both in the UV/optical and the radio/mm 
bands. While curre ntly available measurements have 
moderate sensitivity [111 a nd control of systematics is 
a possible concern—see for a recent overview of the 

theoretical context and observational status—this sensi¬ 
tivity is improving and the prospects are much 

better for the forthcoming generation of observational fa¬ 
cilities fisj] . 

In 0 we extended previously available Principal Com¬ 
ponent Analysis (PCA, see e.g. [l^) methods that had 
been used in other contexts in cosmology to the case of 
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tests of the stability of fundamental couplings, and these 
were then calibrated with some current data from the 
UVES spectrograph in Here we further quantify the 
potential of this method, by discussing the improvements 
in dark energy characterization that can be obtained by 
combining these spectroscopic measurements with those 
of future supernova surveys at various redshifts. We con¬ 
sider two different fiducial models and several a and su¬ 
pernova samples (with different sensitivities and redshift 
ranges) in order to quantify the effect of these param¬ 
eters. In the next section we review the relevant PCA 
methodology, and in Sect. Ill we spell out our assump¬ 
tions on models and data. Sect. IV contains our main 
results, quantitatively expressed in terms of a figure of 
merit, but we also provide a visual illustration by show¬ 
ing examples of dark energy equation of state reconstruc¬ 
tions. Finally we present our conclusions in Sect. V. 


II. PCA TOOLS 

Our study is based on non-parametric PCA techniques 
described in detail in Here we provide a short 

summary of the key features for our analysis. We divide 
the relevant redshift range into JV bins, with the equation 
of state taking the value Wi at bin i, 

N 

w{z) = '^Wi9i{z). (1) 

i=l 

Another way of saying this is that w{z) is expanded in 
the basis 9i, with 9i = (1,0,0,...), 02 = (0,1,0,...), etc. 
The precision on the measurement of Wi can be inferred 
from the Fisher matrix of the parameters Wi, specifically 
from \/(F-i)ji. One can find a basis in which all the 
parameters are uncorrelated by diagonalizing the Fisher 
matrix such that F = KW , where A is diagonal and 
the rows of W are the eigenvectors ei(z) or the principal 
components. These define the new basis in which the 
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new coefficients are uncorrelated, and we can write 

N 

w(z) ='^aiei{z). (2) 

i=l 


The diagonal elements of A are the eigenvalues (or¬ 
dered from largest to smallest) and define the variance of 
the new parameters, a‘^{ai) = 1/Xi. 

We will consider quintessence-type models where the 
same scalar field yields the dark energy and (through a 
coupling to the electromagnetic sector) a variation of the 
fine-structure constant a [11 [H. We take this coupling 
to be 

(3) 

where the gauge kinetic function is linear, 


Brict)) = 1- - (4) 

= SttG and C, is the coupling constant, which in what 
follows will be marginalized over. This can be seen as 
the first term of a Taylor expansion, and should be a 
good approximation if the field is slowly varying at low 
redshift. Then the evolution of a is given by 


/ \ _ ^('^) ^0 ^ f !/ \ ± \ { 

- {z) = -= Ck((/'( 2) - 0o) ■ (5) 

a ao 


For a flat Friedmann-Robertson-Walker Universe with a 
canonical scalar field, (jy^ = [1 -\- w{z)\p^, and therefore, 
for a given dependence of the equation of state parame¬ 
ter w{z) with redshift, the scalar field evolution can be 
written 


(j)[z) - ( 1)0 = 




dz 

1 + z' 


( 6 ) 

where we have chosen the positive root of the solution, 
and therefore the relative variation of a is given by 


— (z) = C r v^3[l + u;(z)]U^(z)^ , (7) 

a Jo ) + z 


where = P(p/ipm + P<p) is the fraction of the universe’s 
energy in the scalar field. 

From the above one can calculate the Fisher matrix 
using standard techniques, as discussed in [l^ for Type 
la supernova data and in Q for fine-structure constant 
measurements. This method has been numerically imple¬ 
mented in a code originally developed by Nelson Nunes 
for 0, subsequently used in , and now modified and 
extended for the present work. We can now proceed to 
discussing our assumptions regarding the fiducial models 
and simulated datasets. 


III. MODELS AND DATA 

We will assume a flat universe, and further simplify 
the analysis by fixing flm = 0.3. This is a standard 


procedure, followed in many previous works including 
SiiSO. This specific choice of has a negligible 
effect on the main results of our analysis. In any case, 
we note that the purpose of these PCA techniques is to 
characterize gains in sensitivity, rather than provide hard 
numbers. 

We will consider two fiducial models for w{z), already 
used in d, [l^ 

wciz) = -0.9, (8) 


i(;s(z) = —0.5-I-0.5 tanh (z — 1.5) . (9) 

At a phenomenological level, these describe two opposite 
and qualitatively different scenarios: an equation of state 
that remains close to a cosmological constant throughout 
the probed redshift range, and one that evolves towards a 
matter-like behavior by the highest redshifts probed. In 
what follows we will refer to these cases as the Constant 
and Step. For each fiducial model we choose a prior for 
the coupling such that it leads to a few parts-per-million 
variation of a at redshift z ^ 4, consistently with m- 
For Type la supernovas we will consider the following 
datasets: 


• A low-redshift sample, henceforth denoted LOW, 
of 3000 supernovas uniformly distributed in the red¬ 
shift range 0 < z < 1.7, with an uncertainty on the 
magnitude of am = 0.11. These numbers are typ¬ 
ical of ’SNAP-like’ future supernova datasets and 
were also used in [l^ (thus providing a useful point 
of comparison). 


• An intermediate redshift sample, henceforth de¬ 
noted MID, of 1700 supernovas uniformly dis¬ 
tributed in the redshift range 0.75 < z < 1.5 and 
with the same am as before. This is representative 
of recent proposals such as DESIRE [l3|, and fore¬ 
casts with a similar sample were already discussed 
(in different contexts) in Bill. 


• A high-redshift sample of sup ernovas identified by 
JWST NIRcam imaging and then character¬ 
ized by extremely large telescopes on the ground 
such as the E-ELT and the TMT. Based on their 
respective Phase A studies we assume a sample of 
50 supernovas in the range 1 < z < 5 for the E-ELT 


2l| and a sample of 250 supernovas in the range 
1 < z < 3 for the TMT (2^ . These will be denoted 
ELT and TMT respectively, and they will provide 
a useful proxy for studying the importance of the 
redshift lever arm versus the size of the sample. 
The redshift distribution of these supernovas is not 
easy to extrapolate, since even the most detailed 
current studies such as those of the SNLS B only 
reach z ^ 1, but in the absence of more detailed 
information we again assume a uniform distribu¬ 
tion in the respective redshift ranges and the same 
am as before. Forecasts with a similar sample were 
already discussed (in different contexts) in [B|- 
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For the fine-structure constant measurements we will 
focus on the ESPRESSO spectrograph for the VLT 
and its successor ELT-HIRES for the E-ELT [ 2 ^; in the 
tables that follow we will denote them ESP and HRS re¬ 
spectively. Both spectrographs include tests of the stabil¬ 
ity of fundamental constants among their key science and 
design drivers. A more detailed discussion and roadmap 
for these tests can be found in [T^. Specifically we will 
consider the following two scenarios 

• A Baseline scenario, in which we will assume mea¬ 
surements in 30 systems with uncertainty CTAa/a = 
6 X 10“^ for ESPRESSO and 100 systems with 
<7Aa/a = 1 X 10“^ for ELT-HIRES, uniformly dis¬ 
tributed in the redshift range 0.5 < 2 : < 4. This 
is meant to represent what we can confidently ex¬ 
pect to achieve from each spectrograph (e.g., from 
Guaranteed Time Observations), given the their ex¬ 
pected sensitivity, and it will therefore provide the 
basis for most of our discussion. 

• An Ideal scenario, in which we will assume 100 
systems with CAa/a = 2 x 10“”^ for ESPRESSO 
and 150 systems with cTAa/a = 3 x 10“® for ELT- 
HIRES. This is optimistic both in the uncertainty 
of individual measurements and in the number 
of measurements. Although several hundred ab¬ 
sorbers are already known where these measure¬ 
ments can be carried out, the sources are quite faint 
and—as can be extrapolated from current VLT 
data 0 —putting together such a dataset would at 
the very least require a very long time and almost 
certainly a dedicated program. Our goal in consid¬ 
ering this ideal case is to obtain an indication for 
the dependence of our results on the uncertainty 
and number of the measurements. 

We note that these are two of the cases already studied 
in 0, whose analysis we will extend. 

For our PCA analysis we will in general assume 30 
redshift bins in the range 0 < 2 : < 4. We note that 
this redshift range is conservative for a measurements, 
since measurements beyond z = A already exist, and the 
infrared sensitivity of ELT-HIRES may push this limit 
even further [^ . In the case where ELT supernovas are 
used, the last bin is extended until 2 r = 5. 

For comparison purposes we will also briefly consider 
a case with only 20 bins. This serves to provide an il¬ 
lustration of the effects of redshift resolution on the re¬ 
construction. Using too few bins is likely to erase useful 
information (especially if the behavior of the dark energy 
equation of state is non-trivial), while using too many will 
lead to very little observational information in some (or 
all) bins, with runs the risk of misinterpreting noise as 
non-trivial information. To a large extent the choice of 
the optimal number of bins will depend on the available 
datasets themselves; while we won’t address this issue ex¬ 
plicitly here, we nevertheless quote the results with the 
two choices of bins to provide the reader with an illus¬ 
tration of the importance of an adequate choice. 


In order to quantify gains in sensitivity we will use two 
different diagnostics. The simplest one is the number 
of PCA modes with uncertainties below apcA = 0.3; 
although this choice of threshold is somewhat arbitrary, it 
has already been used in [ 2 ^ and in our previous analysis 
0 - 

Tables U and |lT] display these numbers, respectively 
for the Constant and Step hducial models, assuming 
the Baseline scenario for a measurements, and for var¬ 
ious combinations of supernova datasets. Considering 
supernova data only, the MID sample adds one mode, 
and the further inclusion of ELT supernovas may add 
an additional one. When combining supernova and a 
measurements, ESPRESSO may add up to two modes 
while ELT-HIRES adds many more. For the combined 
datasets, whether the TMT sample (more supernovas at 
lower redshifts) or the E-ELT one (fewer supernovas at 
higher median redshift) is the more informative one is 
model-dependent. The reason why more modes are well 
characterized in the Constant than in the Step case is 
that, with our choice of high-redshift normalization for 
the a variations, a uniform distribution of the measure¬ 
ments in redshift turns out to be an optimal observational 
strategy for the Constant case, but is far from optimal for 
the Step case; this is further discussed in M- Here we 
have chosen to keep the assumption of uniform redshift 
sampling precisely to highlight this mo del-dependence. 

A somewhat more informative diagnostic is ’figure of 
merit’ defined as the inverse of the product of the un¬ 
certainties of the two best determined modes, FoM = 
l/(o’itT 2 ) [131. We will adopt this for a more thorough 
exploration of these scenarios in the following section. 


TABLE I. Number of PCA modes with uncertainties below 
o'PCA = 0.3, assuming the ’Constant’ fiducial model, the 
’Baseline’ scenario for a measurements, and 30 redshift bins. 



Sne only 

Sne -t ESP 

Sn -h HRS 

LOW 

3 

5 

17 

LOW -I MID 

4 

5 

18 

LOW -I ELT 

4 

5 

16 

LOW -I MID -h ELT 

4 

6 

16 

LOW -I TMT 

4 

5 

17 

LOW -I MID -h TMT 

4 

5 

18 


IV. COMPARING FIGURES OF MERIT 

The figures of merit for the baseline case of a mea¬ 
surements are shown in Tables nil and IlYl respectively 
for the Constant and Step fiducial models. In both cases 
we compare the results obtained with 20 or 30 bins. Ta¬ 
ble |V] shows the results for the Ideal a datasets and 30 
redshift bins, comparing the results for the two fiducial 
models. 

We note that the gains in sensitivity to the dark energy 
equation of state due to ESPRESSO measurements are 












4 


TABLE 11. Number of PCA modes with uncertainties below 
i^pcA = 0.3, assuming the ’Step’ fiducial model, the ’Baseline’ 
scenario for a measurements, and 30 redshift bins. 



Sne only 

Sne -b ESP 

Sne -b HRS 

LOW 

3 

4 

9 

LOW + MID 

4 

4 

10 

LOW + ELT 

4 

4 

10 

LOW + MID + ELT 

5 

5 

10 

LOW + TMT 

4 

4 

9 

LOW + MID + TMT 

5 

5 

11 


relatively modest in the Baseline case, but significant (up 
to about a factor of 2) in the Ideal case. ELT-HIRES, 
on the other hand, will lead to dramatic improvements 
(sometimes more than a factor of 50). These results are 
consistent with the findings of [l^, whose analysis used 
different criteria. 

It’s also noteworthy that judging by this figure of merit 
diagnostic the impact of the E-ELT supernovas is always 
greater than that of the TMT supernovas. This is the 
case whether one is using supernova data only or a com¬ 
bination of supernovas and a measurements. Note that 
for the case of supernovas only, the 50 ELT supernovas 
(uniformly distributed in the range 1 < z < 5) would not 
only be more constraining than the 250 TMT supernovas 
(in the range 1 < z < 3) but also more constraining than 
the 1700 MID supernovas (in the range 0.75 < z < 1.5). 

We caution the reader that there are obvious caveats 
to this comparison. Firstly, finding very high redshift 
supernovas will be difficult, and current estimates of ex¬ 
pected rates are at best uncertain. Moreover, the E-ELT 
or TMT time required to characterize them will certainly 
be costly, and we are taking for granted a temporal over¬ 
lap between JWST and the relevant E-ELT and TMT 
instruments. Nevertheless, the results of this comparison 
do highlight the importance of the redshift lever arm in 
characterizing dynamical dark energy. 

Finally, we also briefly studied how these results will 
be affected if the intrinsic dispersion of the supernova 
magnitudes is am — 0.15 (instead of the am = 0.11 that 
is assumed in the rest of the paper.) In order to quan¬ 
tify this we have repeated the calculation of the figures 
of merit listed in Table m with this degraded disper¬ 
sion, and the revised results are listed in Table ED We 
see that the absolute numbers change significantly (by 
up to a factor of 2 with supernova data alone, by 25 to 
35 percent for ELT-HIRES), but the relative variations 
between the various cases maintain the previously dis¬ 
cussed behavior. Nevertheless it’s noteworthy that the 
figures of merit decrease more for the case of 30 bins, as 
was to be expected. 

One can reconstruct w{z) by keeping only the most ac¬ 
curately determined modes, as was first done, for super¬ 
novas only, in [l^. (See [l^ for an alternative approach.) 
To do this, we need to decide how many components to 
keep. The canonical choice for the optimal value of modes 


TABLE III. Figures of merit for the dark energy equation of 
state, assuming the ’Constant’ fiducial model and the ’Base¬ 
line’ scenario for a measurements, and 30 redshift bins. For 
each pair of entries the top and bottom lines respectively as¬ 
sume 20 and 30 redshift bins. 



Sne only 

Sne -b ESP 

Sne -b HRS 

LOW (20) 

539 

546 

5215 

LOW (30) 

409 

412 

3574 

LOW -b MID (20) 

1090 

1096 

5331 

LOW -b MID (30) 

839 

843 

3655 

LOW -b ELT (20) 

1194 

1215 

8055 

LOW -b ELT (30) 

881 

888 

4947 

LOW -b MID -b ELT (20) 

2371 

2392 

8493 

LOW -b MID -b ELT (30) 

1973 

1980 

5286 

LOW -b TMT (20) 

808 

814 

5302 

LOW -b TMT (30) 

631 

634 

3642 

LOW -b MID -b TMT (20) 

1581 

1586 

5520 

LOW -b MID -b TMT (30) 

1253 

1256 

3814 


TABLE IV. Figures of merit for the dark energy equation of 
state, assuming the ’Step’ fiducial model and the ’Baseline’ 
scenario for a measurements, and 30 redshift bins. For each 
pair of entries the top and bottom lines respectively assume 
20 and 30 redshift bins. 



Sne only 

Sne -b ESP 

Sne -b HRS 

LOW (20) 

536 

541 

1358 

LOW (30) 

404 

407 

982 

LOW -b MID (20) 

1084 

1089 

2003 

LOW -b MID (30) 

831 

834 

1462 

LOW -b ELT (20) 

1225 

1243 

3206 

LOW -b ELT (30) 

881 

885 

1738 

LOW -b MID -b ELT (20) 

2432 

2450 

2561 

LOW -b MID -b ELT (30) 

2175 

2176 

2356 

LOW -b TMT (20) 

821 

824 

1453 

LOW -b TMT (30) 

634 

636 

1055 

LOW -b MID -b TMT (20) 

1605 

1608 

2209 

LOW -b MID -b TMT (30) 

1260 

1262 

1636 


M to be kept corresponds to the value that minimizes the 
risk, defined as [l^ 

N N M 

R = Bias^+Var = ^ {w[zi) - u;*(z,)) Y. 

i—1 i=l j—1 

( 10 ) 

where the notation w means that the sum in Eq. m runs 
from 1 to M and tc* denotes the relevant fiducial model. 
The bias measures how much the reconstructed equation 
of state differs from the true one by neglecting the noisy 
modes, and typically decreases as we increase M. Con¬ 
versely the variance of w{z) increases as we increase M, 
since we will be including modes that are less accurately 
determined. 

An alternative, previously used in a similar context in 
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TABLE V. Figures of merit for the dark energy equation of 
state, assnming the ’Ideal’ scenario for a measurements and 
30 redshift bins. For each pair of entries the top and bottom 
lines respectively correspond to the Constant and Step hducial 
models. 



Sue only 

Sne -k ESP 

Sne -k HRS 

LOW (c) 

409 

996 

58684 

LOW (s) 

404 

554 

11228 

LOW -k MID (c) 

839 

1352 

58737 

LOW -k MID (s) 

831 

955 

11295 

LOW -k ELT (c) 

881 

1515 

79431 

LOW -k ELT (s) 

881 

1064 

18176 

LOW -k MID -k ELT (c) 

1973 

2357 

79639 

LOW -k MID -k ELT (s) 

1971 

2133 

18652 

LOW -k TMT (c) 

631 

1089 

58740 

LOW -k TMT (s) 

634 

712 

11335 

LOW -k MID -k TMT (c) 

1253 

1443 

58846 

LOW -k MID -k TMT (s) 

1260 

1328 

11514 


TABLE VI. Same as Table HVl but now assuming an intrinsic 
dispersion of the supernova magnitudes of (Jm = 0.15 rather 
than (Jm = 0.11 



Sne only 

Sne -k ESP 

Sne -k HRS 

LOW (20) 

303 

309 

1017 

LOW (30) 

229 

232 

730 

LOW -k MID (20) 

613 

618 

1459 

LOW -k MID (30) 

470 

473 

1065 

LOW -k ELT (20) 

768 

769 

1279 

LOW -k ELT (30) 

570 

571 

936 

LOW -k MID -k ELT (20) 

1605 

1606 

1933 

LOW -k MID -k ELT (30) 

1230 

1231 

1450 

LOW -k TMT (20) 

464 

468 

1125 

LOW -k TMT (30) 

359 

361 

807 

LOW -k MID -k TMT (20) 

908 

911 

1541 

LOW -k MID -k TMT (30) 

713 

715 

1130 


is to choose the largest value for which the error is 
below unity, or equivalently, the RMS fluctuations of the 
equation of state parameter in such a mode are 


= ( 11 ) 

and then normalize the errors as suggested in such 
that cr^ = 1 for the worse determined mode and 




1 + a^{ai) 


( 12 ) 


As has been discussed in @ , one expects that the error 
normalization method will generically select more modes 
than risk minimization. This will therefore lead to a more 
accurate reconstruction (ie, closer to the correct fiducial 
model), though correspondingly also one with larger er¬ 
ror bars. The normalization method is therefore the more 


conservative one, while the risk method is more aggres¬ 
sive. We will further quantify this below. 

Importantly, because we truncate the above sum (ne¬ 
glecting the poorly determined modes with high ampli¬ 
tudes at larger redshift) the reconstructed equation of 
state necessarily tends to zero for sufficiently large red¬ 
shift [l^. This unavoidable feature of the PCA trunca¬ 
tion method can be confused with a real increase in the 
equation of state at high redshift. This is another reason 
for wanting to extend the redshift range and the sensitiv¬ 
ity of the measurements: they will lead to more reliable 
reconstructions at higher redshifts. 

Figure [T] shows example of reconstructions, for vari¬ 
ous choices of datasets, for the Step model, using either 
method to truncate the series; these correspond, respec¬ 
tively, to the top and bottom set of 9 panels. In both 
cases we assumed the Baseline scenario for a measure¬ 
ments. Note that in this and subsequent plots we al¬ 
ways plot the reconstructed equation of state until red¬ 
shift z = 4, even though for the E-ELT supernova case 
there will be measurements until z = 5. (In a given figure 
the same mocks have been used for all panels, but differ¬ 
ent mocks were used for different figures; this choice of 
mocks does not significantly affect our results.) 

Naturally the step model provides a best-case scenario, 
since the correct fiducial equation does approach zero at 
high redshifts. This comparison makes it clear that the 
risk minimization method leads to a reconstruction with 
nominally very small error bars, but also one that may 
be somewhat biased. It is clear from the plots that this 
bias can be decreased either by extending the redshift 
lever arm (compare the left-most column, for reconstruc¬ 
tions with supernovas only) or by increasing their sensi¬ 
tivity (compare the middle and right-most set of plots, 
which respectively include ESPRESSO and ELT-HIRES 
a measurements). On the other hand the normalization 
method in this case yields very conservative error bars 
and an almost perfect reconstruction, as long as one has 
a sufficient amount of data deep in the matter era. 

Conversely the case of the Constant equation of state 
provides a worst-case scenario, as shown in Fig. [51 
Here the reconstruction will necessarily be biased at suf¬ 
ficiently high redshifts. In other words, the reliability 
of the reconstruction at high redshift will be model- 
dependent. This much is of course expected, but these 
plots again make it clear that the reliability can al¬ 
ways be increased (all else being equal) by extending 
the redshift range where data is available and improv¬ 
ing their sensitivity—compare the reconstructions using 
ESPRESSO and ELT-HIRES in this figure. Since one 
does not a priori know the high-redshift behavior of the 
equation of state, going as deep in redshift as possible is 
a mandatory aspect of any observational strategy opti¬ 
mization. Tests of the stability of the fine-structure con¬ 
stant can therefore play a key role in this endeavor. It 
remains true in this case that the normalization method 
to truncate the series leads to a more robust reconstruc¬ 
tion. 
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FIG. 1. Examples of dark energy equation of state reconstructions for the Step fiducial model, using the risk minimization 
and normalization methods (top and bottom set of plots, respectively. All panels show the equation of state w(z) plotted as a 
fnnction of redshift. In each set of plots the left panels correspond to supernova data only, the middle ones to the combined 
supernova and ESPRESSO data, and the right ones to the combined supernova and ELT-HIRES data, while the three lines 
correspond to LOW+MID, LOW+MID+ELT and LOW+MID+TMT. In all cases we assumed the Baseline scenario for a 
measurements. 
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FIG. 2. Examples of dark energy equation of state reconstructions for the Constant fiducial model, using the risk minimization 
and normalization methods (top and bottom set of plots, respectively. All panels show the equation of state w{z) plotted as a 
function of redshift. In each set of plots the left panels correspond to supernova data only, the middle ones to the combined 
supernova and ESPRESSO data, and the right ones to the combined supernova and ELT-HIRES data, while the three lines 
correspond to LOW+MID, LOW+MID+ELT and LOW+MID+TMT. In all cases we assumed the Baseline scenario for a 
measurements. 
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FIG. 3. Further examples of dark energy equation of state reconstructions for the Constant fiducial model. Parameter choices 
are identical to those of Fig. [2l except that we now assumed the Ideal scenario for a measurements. 
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In order to further quantify the impact of some of our 
assumptions on the above results, Fig. [3] shows the re¬ 
sults of two other reconstructions for the constant fidu¬ 
cial model case, where we now assumed the Ideal scenario 
for a measurements. This comparison (which obviously 
will only affect the cases with a measurements) clearly 
shows that with the increased sensitivity and number of 
measurements of the Ideal scenario the reliability of the 
reconstruction is significantly improved, both in terms of 
error bars in the various redshift bins and in terms of the 
maximum value of the redshift where the reconstruction 
is not significantly biased. The panels corresponding to 
the Constant model are particularly illuminating in illus¬ 
trating the relative contributions of the supernova and 
a datasets: the difference between the E-ELT and TMT 
supernovas are clearly visible, and indeed even enhanced 
in this case by the presence of similarly sensitive E-ELT 
a measurements. 


V. CONCLUSIONS 

We have used previously available PCA-based forecast 
techniques to quantify the gains in sensitivity expected 
from constraints on the behavior of dark energy enable 
by forthcoming ground and space-based astronomical fa¬ 
cilities. Specifically we have focused on the reconstruc¬ 
tion of the dark energy equation of state, using both fu¬ 
ture space-based supernova surveys in combination with 
high-resolution spectroscopic measurements of the fine- 
structure constant expected from ESPRESSO (whose 
commissioning at the VLT will be in late 2016) and the 
high-resolution ultra-stable spectrograph planned for the 
E-ELT (and currently known as ELT-HIRES). 

Our results quantitatively confirm that the combina¬ 
tion of these two types of measurements, which probe dif¬ 
ferent (but overlapping) redshift ranges leads to a more 
complete—and robust—mapping of the evolution of the 
equation of state, and that a detailed reconstruction be¬ 
tween redshift zero and four is within the reach of forth¬ 


coming facilities. The combination of the two datasets 
leads to figure of merit improvements that are typically 
a factor of a few, and more than 50 in ideal circumstances. 
We also provided a further comparison (for our fiducial 
models) of the two PCA truncation criteria. We would 
argue that the more conservative normalization should 
be the default method for truncation, unless there is 
some prior evidence suggesting that the dark equation 
of state is diverging from -1 at high redshifts (in other 
words, broadly speaking, that one is dealing with a freez¬ 
ing model rather than a thawing one). 

Finally, we should point out an additional possibility 
not directly addressed in our work. A comparison of the 
two reconstructions (separately using supernovas and a 
measurements) can provide a consistency test, and more 
specifically it will be a test for the assumption that the 
same dynamical degree of freedom is responsible for the 
dark energy and the a variation—in other words, that 
one is dealing with a Class I model, in the terminology 
of M- Assuming this is the case, then the coupling 
between the scalar field and the electromagnetic sector 
should be describable (at least to first approximation) 
by Eq. m and one can also infer the posterior likelihood 
for the coupling parameter A detailed study of this 
interesting possibility is left for forthcoming work. 
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